Quasi-uniform in plane magnetization state of thin cylindrical dots in a square array 

and related anisotropy. 



O 

o 
< 

in 



I 

O 

o 



> 
o 

O 



o 

-I— > 



I 

O 

o 



X 
S3 



Konstantin L. Metlov 
Institute of Physics ASCR, Prague, Czech Republic 18221; 
Donetsk Institute of Physics and Technology NAS, Donetsk, Ukraine 8311^ 
(Dated: February 2, 2008) 

The energy (magnetostatic, exchange and Zeeman terms) of a square array of cyhndrical sub- 
micron dots made of soft ferromagnetic material is calculated analytically and minimized, taking 
into account quasi-uniformity of dots magnetization. The dependence of the equilibrium energy of 
the array on the direction of the externally applied magnetic field in the array plane is recovered, 
exhibiting the four-fold anisotropy. The anisotropy constant is calculated. Its values for different 
array geometries are in excellent agreement to the recent independent experiments. New eight-fold 
anisotropy effect is predicted. Theory involves no adjustable parameters. 

PACS numbers: 75.60.Ch, 75.70.Kw, 85.70.Kh 



Advances in fabrication of sub-micron ferromagnetic 
elements and their arrays produced a wealth of experi- 
ments, many of which are still waiting to be explained by 
theory. In particular, the first measurement of the four- 
fold anisotropy of the spin- wave frequencies in square ar- 
rays of circular permalloy dots was reported in 1997ii 

Discovery of the four-fold anisotropy of properties of 
a square lattice might seem obvious at first. However, a 
deeper look reveals that, while the distribution of stray 
fields in a square array of circular magnetic dots de- 
pends on orientation of their uniform magnetization in 
the array plane;^ the magnetostatic energy of the ar- 
ray is completely isotropic for any in-plane magnetiza- 
tion direction.^ It directly follows from Maxwell equa- 
tions that interaction of parallel dipoles (which a uniform 
magnetization state is) depends only on second powers of 
direction cosines and can only produce biaxial anisotropy 
(uniaxial in the array plane) , even if the array is rectan- 
gular. But, as soon as both periods of the rectangular 
array are equal, the (in-plane) anisotropy of the mag- 
netic energy of the uniformly magnetized array vanishes. 

Thus, it is clear that the effect must come from some 
kind of non-uniform magnetization distribution within 
the dots. For example, the two-domain dots (divided by 
the oppositely-magnetized domains in half) do show the 
four-fold anisotropyi^ It can be doubted, however, that 
such a configuration would survive the high (almost satu- 
rating) external field of magnetic resonance experiments. 

Based on the approximate analytical approach to mi- 
cromagnetics of thin flat sub-micron soft ferromagnetic 
elementSjiS the ansatz for quasi-uniform magnetization 
distribution in a circular dot was proposed recently.^ Us- 
ing this distribution as a starting point, the equilibrium 
(in the Ritz sense) energy of the square array of the dots 
in such a state is calculated below, comparing its angular 
dependence to the experiment. 

The magnetization distribution^ can be expressed in 
terms of the complex function of complex variable: 



where z = X/R + lY/R, X, Y are Cartesian coordinates 
on the dot's face, R is the dot radius, i — V—1, Hne 
over a variable denotes the complex conjugation. Com- 
ponents of magnetization unit vector rn are expressed 
through w{z,'z) as mx + imy = 2w/{l + wW) and 
mz = {l — ww)/{l + wW) = 0. The dimensionless param- 
eter p > 1 describes displacement of skyrmions from the 
cylinder's side, their centers (zeros of w{z,^)) arc located 
a,t X — ±pR, Y = 0. The value of p is to be found from 
minimization of the total energy of the array. For p ^ oo 
the magnetization distribution is uniform, for finite p the 
corresponding "leaf" quasi-uniform magnetization state 
is shown in Fig. ^ Compared to the original ansatz in 
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FIG. 1: The quasi-uniform magnetization distribution (0 
with p = 2. Lines inside the particle (its contour is outlined 
by the circle) are tangential to the magnetization vectors. 

Ref. 6 the expression ^ allows the average magnetiza- 
tion of the dot to point in arbitrary direction at an angle 
ipo w.r.t one of the array axis (this was unnecessary com- 
plication for a single doli^). 

To evaluate the magnetostatic energy of the array of 
such dots the magnetic charges formalism is used. It 
is based on the assumption that the problem is static 
and there are no currents in the system. This allows 
to express the stray magnetic field through a scalar po- 
tential, created by a system of "fictious" (introduced for 
mathematical convenience) magnetic charges. The vol- 
ume density of these charges is av = — div(m). On the 
boundary of magnetic material cry reduces to the surface 
charge with the density as = {rn ■ n) proportional to the 
(boundary) normal component of the magnetization vec- 
tor. Both densities are normalized here by the saturation 
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magnetization of the dot's material {Ms). Inside of a dot 
with magnetization distribution ^ they are: 



av = 



2r cos{ip — ipo) 



^r** + p"* — 2p27-2 cos 2(<y9 — ipo) ' 

(p^ - 1)C0S((^ - If a) 

^l+p4_2p2cos2(^-(^o)' 



(2) 
(3) 



where r < 1 and < < 27r are the polar coordinates 
on the dot's circular face, normalized by the dot radius. 
These are the rotated charge densities of Ref. 

To fully take into account periodicity of the array and 
all the implied interactions between dots, the density of 
magnetic charges is represented by its Fourier compo- 
nents: 




ipA sin fc^izTT 
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F = 



^ — 27Tip{nX-\-rnY) ^ — 2-Kipr{n cos ip-^rns\i\i^) 



(5) 
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where p = RjT , X ^ Lz /T , T > 2R tmd Lz are the 
array period and thickness, ao,o = by construction. 

Because the Fourier basis functions F are also the 
eigenfunctions of the Laplace differential operator, the 
Poisson equation (resulting from the Maxwell equations 
under assumptions of the magnetic charges formalism) 
for the scalar potential becomes algebraic for its Fourier 
harmonics and can be readily solved. Normalized mag- 
netostatic energy per array cell volume due to the in- 
teraction of the stray field with the magnetic charges 
GA/ = EmKp.oIbMIT'^Lz) is then 
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where |a„,,„|2 = a„^„ia_„^_,„, C = ^z^z, A^o is perme- 
ability of vacuum in SI units, 7s is a units-dependent 
factor equal to 1 in SI (in CGS /^o = 1, 7b = 47r). The 
integral can be taken, allowing to represent 
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To evaluate the 2-d Fourier components let us first 
note that the denominator of both surface and volume 
charges can be expanded into the power series in (a 
small parameter in case the magnetization state in the 
dot is quasi-uniform) using the identity for the generating 
function of the Legendre polynomials Pi{x) 
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— 2p^r^ cos tp 
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where ip = 2{(p — (pQ). Using the identities for the gen- 
erating functions of Bessel's functions (formulas 6.521 in 
Ref. 0), the Fourier basis ^ can be represented as 
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where a = 2TTipr and Jk{x) are the Bessel's functions of 
the first kind. 

In such a representation it is easy to evaluate the an- 
gular integral ©, which, for a given i, is non-zero only 
for fc -I- Z = ±(2j + 1) with j = 1, 2, . . . , i. This allows 
to remove (by making it finite) one infinite summation 
from a triple i, k, I sum, resulting from expansions H10|l . 
(|ll|l . Then, using Bessel's summation theorem, each of 
the remaining infinite sums over k + I = (2j -I- 1) and 
k + l = —[2] + 1) can be evaluated and combined, yield- 
ing the representation of a„^m term for a given i in terms 
of finite sum of Bessel's functions of the orders 2j -I- 1 for 
j = 1, 2, 3, . . . , «. The resulting terms can be easily inte- 
grated in r for the part of the integral representing 
the volume charges. This calculation is quite voluminous 
in general.^ Fortunately, to obtain the expansion of the 
magnetostatic energy up to the 1/p'' (the first neglected 
term is ^ 1/p^) it is sufficient to evaluate the volume 
charges Fourier components for i = 0, 1 and the surface 
charges for i = 0,1,2 (because of p2 jn Eq. 13). 

The final result for the magnetostatic energy can be 
represented as 



Cm = eo + 



ei + cos(4(po)eia , e2 + cos(4(^o)e2a 
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Kiix 

K2{X 
K2a(x,tlj 



Jl{x), (14) 
Ji{x)J^{x), (15) 
Ji(a;)J3(2;)cos(4V'), (16) 
-~2Ja{x)Ux)/x, (17) 
( Ja {xf + 2 Ji (x) Js {x)) cos(4V')/2. (18) 



Functions — ei{p,X) depend only on geometry of the 
problem and take into account the self-energy of the vol- 
ume and surface charges, as well as interaction between 
them, both in a single dot and across all the dots in the 
array. 

These expressions can be verified by renormalizing 
them to the unit of magnetic dot volume (as opposed to 
the volume of lattice cell) and taking the limit T oo. 
In this case eia — e2a — 0, and cq, ei and 62 become 
identical to the corresponding functions e", 62 and 64 of 
Ref. 0. Because the same thoroughly verified automation 
procedure was used in deriving all the functions Ki, in- 
cluding the angular ones (marked with "a"), it can be 
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expected with a high degree of confidence that all the 
expressions for ei{p,X) are correct. 

Looking at 1)12(1 it is immediately seen that for uni- 
formly magnetized dots (whenp — cxd) the magnetostatic 
energy shows no angular dependence. The function cq is, 
thus, the energy of a square lattice of uniformly magne- 
tized dots, which can be independently verified. 

Compared to the above, evaluation of the exchange 
and Zeeman energy terms is trivial. Both contain no 
interaction between dots and can be evaluated for each 
dot separately. Actually, both were already calculated 
in Ref. 0. Renormalizing (2) from Ref. as e^; — 
EE/if^olsMgT'^Lz) and using the expression for the av- 
erage magnetization (16) therein we get up to 1/p*: 

(19) 

where t — T/Le, Le = \JC/{p,QMg) is the exchange 
length of dot's material, h = H/ (ho^bMs) is the normal- 
ized field magnitude, and tpi is the angle of the in-plane 
applied field with respect to the lattice. 

Now it remains to minimize the sum of ((1211 and 1(1 9() 
to find the equilibrium values of p and ipo. This results in 
two solutions for (pQ separated by a small angle from the 
direction of the applied field tpi , meaning that there will 
be hysteresis, depending on whether the current direc- 
tion of the applied field was approached from bigger or 
smaller angles. Because the typical magnetic resonance 
measurement is performed in high sub-saturating mag- 
netic field, this hysteresis can be expected to be small. 
Let us neglect it here for simplicity by putting (pQ = (pi. 
Minimizing we get 

p2 ^ feTTjO^ + 12(7r/(7Br^) + 62 + cosAipi) ^^^^ 
6(-ei - eiQC0s4<^i) 

Let us assume for now that the denominator here is al- 
ways positive (that is ei < 0, eia > and jeial < |ei|), 
which is the case for all the experiments analyzed below. 
Breaking of this assumption will be discussed at the end. 
The value of p completely defines magnetic structure of 
dots in the array. Substituting it back, the equilibrium 
energy can also be expressed straightforwardly, but let 
us now concentrate on its angular dependence. 

The largest (besides the constant term) Fourier har- 
monic in the angular dependence of the equilibrium en- 
ergy is proportional to cos 4(^i , which can be clearly clas- 
sified as the four-fold anisotropy. The higher harmonics, 
while modifying slightly the dependence itself, do not al- 
ter positions of extrema on the angular dependence with 
minimums located at (pi = 7r/4, 37r/4, 57r/4, 77r/4 and 
maximums at ipi = 0, 7r/2, tt, and 37r/2. Let us, thus, 
introduce the fourfold anisotropy constant, representing 
the equilibrium energy of a single dot (that is, normal- 
ized by the dot volume, which is the volume of magnetic 
material in the array cell) as E — Eq + K4 cosAipi with 
K4 > 0. The value of K4 is the half of the equilibrium 
total energy (renormalized by the dot volume) difference 



between configurations at = and ipi = 7r/4. Further, 
noting that in the resonance experiments the high exter- 
nal field plays a more notable role in stabilization of the 
"leaf" state than the exchange interaction (stabilizing it 
in the absence of the applied field), a limit of t — > 00 is 
taken to simplify the expression into 

^ _ fio-/BM^ 3{ei-eiaf 3(ei + ei„)^ \ 

(21) 

The full expression, more accurate for the case of the 
array period comparable to the exchange length, can be 
easily obtained by the reader (it is also plotted in figures 
here). Nevertheless, for the subsequent comparison to 
experiments the precision of ((21(1 is sufficient. At large 
applied fields H/{p,qJbMs) » 1 it is possible to expand 
K4 asymptotically as 

A,.-„.™M|5i£| + o(Jj)^ (22) 

Let us now apply the obtained knowledge to describe 
experiments found in literature. The first and foremost 
is the experiment found in Ref.Q. To interpret it, the au- 
thors added the anisotropy of the form sin^ -0 cos^ ip 
(tJj is the angle between the magnetic moment in the dot 
and the lattice axis) into their numerical program for 
finding the resonance frequencies and fitted the angu- 
lar dependence of the spectra to determine the K^. Be- 
cause the magnetization of the dot is close to uniform 
(-0 = (pi) this term produces the shift in the total dot's 
energy, followed by the measured spin-wave mode. Due 
to sin^ (^1 cos^ oscillating between and 1/4, while 
cos4(/3i between —1 and 1 with maximums of one cor- 
responding to minimums of another (and vice versa), 
there is additional factor —8 linking the anisotropy in 
Ref. I' to the definition here. Another peculiarity is that 
the value of in Ref. Q tends to a constant value at 
large applied fields. This can not be explained by any 
model of the quasi-uniform state of infinite array, which 
at /i — !■ 00 must transform into the uniform state, show- 
ing no anisotropy. Thus, it seems another factor (proba- 
bly due to the array shape effect) is present in the exper- 
iment. Assuming this factor is field-independent, let us 
simply add it, expressing — —8K4 — 0.6 • lO^erg/cm^, 
plotted in Fig. [21 (taking Mg = 850 emu, typical for 
permalloy, but not specified in Ref. T) along with the 
experimental data. The agreement seems to be rather 
good, noting that the precision of the experiment did 
not allowi to reliably resolve between dependencies of 
the Kl for two considered arrays. 

A much more precise recent experiment was done us- 
ing FMR technique. 12, In FMR it is impossible to change 
the value of the applied field, which was fixed at llOOOe, 
but the authors tracked the anisotropy field as a func- 
tion of the array geometry, measuring different arrays 
with periods from 1100 to 2500 nm. The anisotropy 
field can be expressed through the anisotropy constant 
as H4 = 2K4/MS, which is plotted in Fig. 13 along with 
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r=1100nm ^ - 

p=(1000 nm/2)/r=0.45 
Mj,=850emu, L^=20nm- 



2 4 6 8 10 12 14 16 

H/ (471 Mg AO 

FIG. 2: Dependence of the four-fold anisotropy constant K4 
on the magnitude of the in-plane magnetic field. Points are 
from Ref. U The solid lines l|21|l and the full expression 
(without taking r — > cx) limit) are indistinguishable. Dashed 
hues are from Eq. 1221 The upper set of lines corresponds to 
Lz = lOOnm, the other to Lz = 50nm. The values of the 
in-plane demagnetizing factor A'^ are shown in the legend. 
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FIG. 3; Dependence of the four-fold anisotropy field H4, on 
the inter-dot spacing. Points are from Ref. The solid line 
plots the Eg. 1211 (again, the finite r case is indistinguishable). 
Dashed and dotted lines show the asymptotes of II21II up to 
0{l/h) (Eq.|m and 0{l/h^) respectively. The value of Le 
obtained from the measured— exchange constant (A = C/2). 



the array parameters and the experimental data from 
Ref. ^3- The agreement is excellent for the case of K4 
calculated from H21|l . whereas the simplified asymptotic 
expression (|22ll and its higher order equivalent are clearly 
not good enough at such a small field. 

Let us now look back at H20() . It turns out that it 
is quite possible for the denominator to become zero and 
negative. This happens for dot geometries approximately 
defined by the numerical fit Rcr/T = 0.51-0.37-{LjT) + 
0.13 • {L:,/Tf in < L^/T < 1 with fixing the result- 
ing Rcr/T to 0.5 when it exceeds this value. At radii 
above Rcr the parameter p becomes imaginary for some 
directions of the applied field, which corresponds to the 
transition to the so-called "flower" state. In this state (as 
opposed to initially considered "leaf" ) the magnetization 
vectors diverge outwards of the line through the dot cen- 
ter parallel to the applied field. All the expressions here 
still apply in the case of "flower" (with imaginary p), but, 
because of the additional "flower" ^ "leaf" transition the 
dots array with the R > Rcr shows the additional eight- 
fold (!) anisotropy of the form K^cosAipi + ifg cos 8(/3i 
with Ki still given by (|21|l and < 0. This transi- 
tion and the resulting anisotropy will be explored in the 
forthcoming extended paper on the subject. 

To conclude, the considered problem belongs to the 
class, where the leading order energy contribution van- 
ishes due to symmetry, and the system is governed by 
higher order energy terms. In most other cases details of 
the quasi-uniform magnetization distribution make only 
small (and often negligible) corrections to the leading 
order results, but here they are responsible for the com- 
pletely new effect, absent in the leading order. Apart 
from quantitatively describing the four-fold anisotropy 
in the magnetic dot arrays, the presented theory asks for 
experiments on thicker dots, which should reveal the new 
eight-fold anisotropy effect. This theory can be expected 
to describe it quantitatively. 
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